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High energy hard scattering in large Nc hmit can be described by the QCD dipole modeL In this 
paper, single, double and triple BFKL pomeron exchange amplitudes are computed explicitly within 
the dipole model. Based on the calculation, the general formula 7q*'* = (^X (§)) which governs 
the anomalous dimension of 1 fc amplitude is conjectured. As far as the unitarity problem is 
concerned, we find that the anomalous dimension 7 varies from graph to graph due to the DGLAP 
evolution. In the end, a comparison between this computation and reggeon field theory is provided. 

■ PACS numbers: 12.38.Cy; ll.lO.Hi: 11.55.Bq 
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Q I. INTRODUCTION 

:z;: 

High energy QCD small-x evolution of hadron can be described by the Bahtsky-Fadin-Kuraev-Lipatov (BFKL) 

■ Pomeron [l], [3, 3- BFKL pomeron should be relevant for describing the exponential growth of scattering amplitudes 
I with respect to rapidity. More than a decade ago, the dipole picture of pomeron was discovered by Mueller et al]3,[l,[l], 

r~~i- and the exact equivalence [t'I between the color dipole model and the BFKL pomeron result was verified several years 

' later. Onium-onium scattering at high energy is used to describe this dipole picture of high energy hard scattering 

II in large Nc limit. Supposing that M, the mass of the onium, is large enough to justify the fixed coupling calculation, 

CIh one can find that the onium-onium scattering cross section at single pomeron exchange level behaves as 

; 

r^/ a ^ xx'a^exp[{ap -1)Y] (1) 

■ where x, x' are the sizes of the interacting onia, Y — \v^ and ctp — 1 = SqCf in 2 .^j^j^ (j^ _ ~2N~' ^ large A^c 
limit. Besides the exponential growth of the cross section, one also gets the anomalous dimension 7 = ^ as a result 
of a being proportional to x2:'.(The shift of the exponents of x^ and x'"^ from 1 is called the anomalous dimension of 
the BFKL pomeron.) 

T-H On the other hand, single pomeron exchange amplitude violates the unitarity and the Froissart bound at extremely 

, high energy. Multiple exchanges of pomerons (here we mean pomeron loops) should be able to reduce the growth 
' rate of the amplitudes, and eventually unitarize it at high energy. Within the dipole model, Kovchegov equation^, a 
f — ■ nonlinear evolution equation, is derived by re-summing the fan diagrams(i.e., multi p le p omeron exchanges in dipole- 
nucleus scattering). Recently, there has been a lot of development P. [lOl [Til [T^ . [Tsl . fisl [igI in high-energy QCD 
^ \ evolution. Evolution equations which include pomeron loops are derivedjlTI. llSl. llOl. |2C|. l2ll |22| and utilized to QCD 
phenomenology (e.g., see [13, US, H^, H^-)- far as the unitarity problem is concerned, one has to integrate over 
rapidity and stretch the pomeron loops as large as possible in rapidity space in order to obtain the maximum leading 
order loop amplitudes. Meanwhile, the anomalous dimensions of pomeron loops are determined by the lower and 
upper pomerons with infinitesimal rapidity length. Assuming that only one pomeron can interact with the target 
or projectile onia, these two pomerons then connect the target onium and projectile onium, respectively. In spite of 
having the same characteristic function as other large rapidity pomeron, these two pomeron, however, give rise to 
different value of anomalous dimensions to the loop amplitudes [sl. IgI. [27| . 

As building blocks of pomeron loops, 1 =^ 2 amplitude as well as triple pomeron vertex, which features dipole 
pair correlation in onium states, has been studied extensively [1, [13, [1^ [H, [l^i [U, [13 ■ C)ne can easily obtain the 
anomalous dimension of one loop amplitudes by calculating the 1 2 amplitude Uy since the one loop diagram 
is just nothing but square of 1 2 by hooking them in the middle of the rapidity. Following the same reason, the 
1 3 amplitude indicates the anomalous dimension of two loop amplitudes. The objective of this paper is then to 
investigate the 1 =J> 2 and 1 =J> 3 amplitudes explicitly after integrating the rapidity from to the maximum rapidity 
of the system Y , and then generalize the results to 1 fc amplitudes. 
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In this paper, in the way described above, we expUcitly calculate the anomalous dimensions of the single, double 
pomeron and triple pomeron exchange amplitudes in QCD dipole model in the leading logarithmic approximation. 

In addition to the usual pomeron anomalous dimension ^, we find that = 0.18 and 7^^ = 0.12 for double 

and triple pomeron exchange amplitude, respectively. The anomalous dimension 7q = 0.18 of one loop amplitude 
is not new, and it has been obtained in ref. @, [23|. The anomalous dimension of the triple pomeron exchange 
amphtude 7q'^''* = 0.12 is new and is given by Eqs. ([79]) and ([80]) . Based on the calculation, the general formula 
7^,'^^* — {kx (^)) which governs anomalous dimensions of 1 ^ fc amplitude is conjectured in the region In^ (gp) <C 
i4aCFC(3) Y p2q2 ^ -|^^ where p stands for all dipole sizes and q represents all momentum scales. 

The paper is organized as follows: we start with the computation on single pomeron exchange amplitude in onium- 
onium scattering and discussion on the DGLAP evolution in double logarithmic limit, then calculate the double 
and triple pomeron exchange amplitudes, which are then followed by the generalization to 1 A; amplitude and a 
comparison with the reggeon field theory, as well as conclusions. 



II. THE DIPOLE DENSITY IN THE QCD DIPOLE MODEL. 
A. Single BFKL pomeron exchange in onium-onium scattering. 



In order to be intuitive and complete, let us first sketch the well-known single pomeron exchange amplitude in 
onium-onium scattering. Following 0,0,11^, the distribution of dipoles in an onium state in coordinate space reads. 



n— — oo 



2v'^ + f duj 



2iii LU — uj {n, v) 



ra 



where all the coordinates p are complex coordinates in the two dimensional transverse space, 

~P = {Px,Py) 

P = Px+ iPy and p^ p,,- ipy 



(2) 



(3) 
(4) 



E ■ [Po-yPi-y) is eigenfunction of the SL (2, C) group, and 



E' 



Poi 
Pa-fPi-f , 

17P17 ) 



Poi 

POlPl-i 



with p„ = - p^^ , = po - and 



■ 

n = — h 11/ 



2 

1 + n 



IV. 



(5) 
(6) 



(7) 



We define 



w {n, v) 



AaCp 



1 , fl 



(8) 



where ip {x) is the digamma function. For fixed values of one can easily find that lo (0, v) is always larger than u) (n, u) 
with nonzero values of n. Thus, as one can find in later discussions, the n = contribution corresponds to the dominant 
pomeron trajectory with a positive intercept, and 71 ^ parts correspond to sub-dominant reggeon trajectories with 

intercepts being equal to or less than 0. For convenience, we also define x(i^,7i) — '0(1) — ^V' ^ ^^j"^ ^'^) ^ 

(y^-^-^ — iv^ , where x ('^i is a real function of v which is analytic in the strip — ^-^^ < Im v < ^^y^- ^'^^ '"^^l 
value of v, x ('^j is symmetrical with respect to = 0, and it has a global maximum value at = 0. 
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According to the Fourier transform of E'^''^ (P07P17) 



with R=^{po + pi) - p p = Pq- and 



7732^''' r (|n| + 1/2) r (|n| /2 + iv) 

\n\ /2-ivT {\n\ /2 + iv + 1/2) V (|n| /2 - iv) ' 

(1) 



One can get rid of the impact parameter dependence and define ny' in momentum space, 



^ f ^ 2iy^ + n^/2 1 ^ 2,(2), 

/■ c^g; e"^ bn^.K., \ppj ^n,.. ^ 

^7 2^^0,-0. (n,.) (2.2)2 (PJ- 



(9) 



(10) 



(11) 



(12) 



where r = 5 (Po + Pi) ^-nd ra 



- Poj). At very large rapidity Y, the term corresponding to n = dominates 



the exponent since uj{0,i>) is greater than other Lj{n,i>). Thus, all the contributions of sub-dominant trajectories 



(contributions from nonzero n) can be neglected. Therefore, defining ny^^^^q^ 
that 



n 



bo,v \p\ 
27r2 



£^^(P) = / 



-1 l+2iv 



\R--2\ 



47r^(5^^^ (g — ga) , and noting 



(13) 



one gets 



''Y,q 



2v^ 1 bo,„bQ \ppj ^C10,!^*/ ^ t rn 

dv—r — E'^q'" (p) E''^'" (p„) exp [uj (0, z/) Y] 



4 2 
Pa 



{2n^y 



1 kpl , , , , [(ap - 1) ^] cxp 
re (p, g) e (p„, -q) 



^1" (P/Pg) 

28qCfC(3)F 



2 kPal 



(7aC (3) Ci.r) 



3/2 



(14) 



(15) 



with the famous BFKL pomeron intercept ap — 1 = 8«gF in 2 ^ ^j^^j ^j^^j being much smaller than 1. In 
arriving at the above result, we have used saddle point approximation and assumed that diffusion approximation 
In^ (p/p„) «; i4aCFC(3) y jg ^^YiA, and have defined 



1 /■ 1 



R + 



(16) 



It is straightforward to see that the saddle point approximation picks the dominant contribution of the J d,v integral 
in the vicinity of = which gives rise to the anomalous dimension of single BFKL pomeron: 7* = i + iv* — i. 

Furthermore, employing saddle point approximation, it is easy to find a general formula for the intercepts of all 
reggeon trajectories(for all values of n): a„ — 1 = ^"^^ X (0) '^)- For reader's convenience, we list the first three sets of 
intercepts in the following: ao — 1 = 2.77-2^^, a±i — 1 = and a±2 — 1 = — 1.23 ^"'-^^ . The n = contribution gives 
to the dominant pomeron trajectory with a positive intercept, and n 7^ parts correspond to sub-dominant reggeon 
trajectories with intercepts being less than or equal to 0. 

In addition, the single pomeron exchange amplitude between two onia with sizes p and p' then scales as 



(p,p',g,Y)oc(gV)M9V^)" 



^2 "xp [{ap - 1) Y] 
' (7aC(3)C^y)'/2 



(17) 



1 Hereafter, we use p as an abbreviation of above definition for pQ — p^. 
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B. The DGLAP evolution. 



In comparison, we would like to discuss another interesting limit of the one dipole amplitude. Other than the 
extremely large Y rapidity limit, we now focus on the limit where In (^^^ ^ ia^j^y and In ^ ^ 1 while Y is 

relatively small. This calculation is useful to understand the new anomalous dimensions of dipole pair and triplet 
densities found in later discussions. Therefore, one can cast the one dipole amplitude fEa. (fT4)l ) into the form. 



n 



(i)p,p„ 



d^i^i??-- (p) E^* (pj exp [io (0, v) Y] 



Pi 



dv- 



27r2 



9 



exp [lo (0, v) Y] , 



PI. 



iaCp 



Y 



^21n(^) ^ln(^) 
where one has used 6o,,y65 — ^ and the expansion 0, [111 



: exp 



(18) 
(19) 

(20) 



(p)£;°^- (pj 



9<7 



as well as saddle point approximation in the vicinity of 7 ~ 1 or 0. Thus the cross section scales as, 



^■(P.Pa:^) ^oi p exp 



(21) 



(22) 



In this coUinear limit (p^ 3> p), we retrieve the well-known result of the DGLAP evolution (3^ when the strong 

coupling is fixed. The DGLAP evolution is featured by the exponent being 2^/ '^"^^ Y In ^ since it involves 

evolutions both in rapidity and virtuality (the reciprocal of the dipole sizes) in the so-called double leading logarithmic 
limit. Furthermore, as a result of being proportional to p^ when the logarithmic dependence of p^ is neglected, the 
anomalous dimension in this situation is zero, as opposed to the anomalous dimension of the BFKL pomeron being 
^(for a pedagogical introduction of BFKL and DGLAP evolutions, see [35]). In addition, in small Y limit, this 
coincides with the perturbative QCD calculation of the dipole-dipole cross section[3, 



add {x, x) = 2T:a^x\ ^1 -I- 



In- 



x< 



(23) 



with being the greater one of x and x' and a;< being the lesser. Therefore, one expects that the anomalous 
dimension should approach zero when rapidity is small or DGLAP limit (here we mean double leading logarithmic 
limit which includes both evolutions in virtuality and rapidity) is valid. 



III. THE 1 ^ 2 AMPLITUDE AND DOUBLE POMERON EXCHANGE. 

As shown in Fig. [Tl the 1^2 amplitude is defined as the dipole pair density in an onium state. It involves a 
triple pomeron vertex which splits the upper pomeron into two descendent pomerons. In QCD dipole modeljl, @], 
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FIG. 1: Dipole pair density ny in an onium state with initial size p = pg — and two children dipoles having transverse sepa- 
ration and pij. The left graph shows the rapidity structure of n2, and the right graphs indicates the graphical representation 
of the triple pomeron vertex. 



following Peschanski|28l|. one can write the 1^2 amplitude as 

n-Y^ {PoPl, PaoPai^ PboPbi) 



dhdhndhh 



dhJ 



PaPh J 27ri LU — LU [n,i>) LU — LU (a) — uj (6) 



X / d'p^d^Pf.d^p^E'^-'' (po^,Pi^)i?''-''' {Paoa^PajE'^'^'"' {PboP,Pb,p) 



d^p^d'^p^d'^p^ 



IP23P34P42I 



2 (p2,,P3.)£^''-'^*(P2a,P4a)£^ 



■hhMh 



[P3P'P40) 



dhdhndhh 



1 



pIpI Jo 



aNr_ 



dy exp [uj {n, v) {Y ~ y) + [u (a) + uj (b)) y] 



X / d^p^d'ppd'p^E'^^'^ (P07,P1^) i?''"'"' (Paoa,Paia) E'^'^'"' {PboP^Pb.p) 

d^P2d^P3d^P4 j^h.h* 



|P23P34P42r 

with p„ = p^^ - p^^ , pf, ^ pf,^ - and 



IE" 



UJ (a) 

L0{h) 



AaCp 

■K 

AaCp 



■ 1 /l + |n,| 



1 , fl + \na\ 



TT 



1 , A + 



+ ivb] - 2^ 



1 , /l + l"6| 



We also use a compact notation, 



4 ■ 



Essentially, the above expression for is the same as the ones in ref. [29|, |30|. On the other hand, hereafter, we 

(2) 

employ different methods of evaluation and show that there exists a different anomalous dimension in rty other than 
the typical pomeron anomalous dimension 7 = 1/2. 

(2) 

First of all, the last line of the above definition of Uy can be defined and written as 



(24) 



(25) 



(26) 



(27) 



—h,ha,hb 



R. 



R. 



l-hA~ha,l~hb 

d^p^d'^p^d'^p^ ^^jl 



E'^^'^* (p2,P3.) i^""'""* (P2.P4J e'^'"'* (P3/5P4/5) 



IP23P34P42I 
After using the SL (2, C) transformation 



P 



, _ Pf3- P-i P- Pa 
Pf3- Pa P- P-y 



(28) 



(29) 
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one can easily write 



-+i -1+fcA (--~i+l--i+]--i+k 



where 



and 



g3pih,ha,hb 



i — -h + ha + hb, 
j = -ha + hb + h, 
k = —hb + ha + h, 



d^P2d P^d Pa h ( P24 



(30) 



(31) 
(32) 
(33) 



|P23P34P42r 
P34 



P23P23 



P2P4 



P2A 
P2P4 



P34 



(34) 



[l- p^){l- p,)) V(l-P3)(l-P4) 

Explicit evaluation of g^^p can be found in ref. ([aH) and Appendix E 

Furthermore, according to the Fourier transform of E'^'^ (pp^pj^^), one can get rid of the impact parameter depen- 



dence and define riy •* in momentum space 



{PoPl^PaoPai^PboPbi] 



duj 



dhdhadhb^^2 / r, • r \ / \ fi\ 

PaPb J UJ — UJ (n,h') LO — Lo (a) — uj (0) 

bn,vbna,Va^ni,Ui \PPaPb\ T^n.v 



(2^2)o 

X / d^p^d^ppd^p^ exp {iQp^ + iqaPa + ^^ftP/s) 



i?"-"^ (p) K;'"-" (pj K-^^ (p,) 53P {l-h,l-K^~ hb) 



-l+fc 
Paf} Pp-y Pia 



\a.h. 



(35) 



where a.h. stands for the antiholomorphic part of the square-bracket term, r — j {Po + Pi); '''o. — \ {Paa + Pai) ^-^^ 

fb = \ {Pbo + Phi)- Hereafter in this section, we compute riy-* in p^q^ ^ 1 limit, where g^, and are of the same 
order and p, p^ and Pb are of the same order as well. 
In addition, let us define 



Q,qcL,qb 



In order to evaluate Ig'ql'q^, we change the variables p^,, p^ and p^ into w, v and w, where 

^ = ^ (Pa + P/3) ' 



d^p^d^Ppd^p^ exp (iQp^ + zgaPa + ^gbP/j) 



Pap Ppi P-ya 



(36) 



^' = Pa-Pp^ 
W = Pa- P^+ Pp- Pj- 



(37) 



Then the integral can be cast into, 



- / d^ud^vd^w exp [i {Q + qa + qb) u - iQ— + i {qa - qb) 

[a.h.] , 



W + V 



-l+fc 



7r2(5(2) (Q + g„ + g,b) / d^-yd^wexp (~iQ ■ ^ + i {qa - qb) ■ ^) 



W — V 



-1+3 



W + V 



-l+fc 



[a./i. 



(38) 
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where the factor of j comes from the Jacobian. Hereafter, in order to simplify the calculation, we define 

with Q- — ^ [qa — qt), perform the Taylor expansion of the term exp {—iQ^) and only keep the first term. Keeping 
only the first term is equivalent to the physical case when Q = and qa = —qt = q which corresponds to the forward 
scattering of an onium on two nucleons. We put the evaluation of Iq^'^'^'' and discussion of higher order terms in the 
expansion in Appendix |X] and |B1 respectively. As shown in the appendix, the following conclusion holds for all the 
terms in the expansion when n = 0. 



Therefore, one obtains. 



n 



1 r aNr 



YQ nO = / dhdhadh^ / -dy exp [lo {n, v) (F - y) + (c^ (a) + (6)) y\ 

PaPb Jo 



{27T-) 



x4 2 /^Y+^+'=+*+^+'=-^ r(fc)r(j) r(i-J-fc) r(z + j + fc-i) 

\q-J r(fc + j) r(i-J)r(i-fc)r(2-z-J-fc) ^ ' 

where the last subscript of g'^^'g'' means that this result comes from the first term of the expansion. Assuming that 
the diffusion approximation is valid which requires In^ (qp) <^ i4aCFC(3) i^f^ f^j-st evaluate the rapidity-integral. 
The integral yields 

{exp [{u; (a) + lu (5)) Y] - exp [lu (n, iy)Y]} . (41) 



u) (n, u) — u) (a) — uj (b) 



For the first term, the saddle point approximation fixes ha = hb = na = nt = 0- As a result, ny is then proportional 

(2) 

to exp [2 {ap — l)Y]. In comparison, for the second term, Uy' is proportional to exp [{ap — 1) Y] because saddle point 
approximation fixes h = ^ and n = 0. Therefore, we can drop the second term in later discussions since it is the 
next-leading order contribution which is exponentially suppressed. 

In the following, we break the discussions into two parts. The first part is on the rt = case, and the second 
part discusses the result when n ^ 0. Usually, we can neglect the contribution from sub-dominant trajectories of 
pomerons when rapidity is large enough as we explained in the calculation of the single pomeron exchange amplitude. 
Nevertheless, the rapidity interval of the upper pomeron after integration is infinitesimal in this situation. Therefore, 
we have to seriously consider the contribution from sub-dominant trajectories. 



A. n = part 



n = part of the amplitude is the only angular- independent part. From Appendix [Cl one uses saddle point 
approximation to evaluate / dhadhi, integrals. 



dh„ 



^ria.Va \Pa\ J^na,Va 



27r2 



iPa) 



q- 



I p exp [{ap - 1) Y] exp 

^7' ^'"''"^ {7aa3)CpY) 



exp [oj (a) Y] 

28qCfC(3)V 



3/2 



(42) 
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FIG. 2: Illustration of the contours in complex 7 plane for the computation of riy when n = 0. 



where 7a = 5 (^a + ha) 



i + Wa- In the end, there is only / dh left in the expression, 



Eq {p)93P ( 1-7, 2' 2 



dh 
x47r^ 

1 PaPb 

Let us focus on the dft,-integral (first two lines of the above expression) 
47r2 f „ 1 



(7)r(i-27) 
r(27)r3(i-7) 

exp [2 (ap — 1) y] exp 

e-{Pai<la)e{pf,,qb) ■ 



^(ln''(gaPa)+ln''('?t.Pt)) 
28qCfC(3)V 



(7aC (3) CfYT 



bo,u \p\ j^o,u I (. 11 



dh- 

d'y 



4 r3 (7)r(i-27) 



r(27)r3(i-7)' 



2 /• d7 IpI . A 1 IV 4 



7-1 



r(7) 
r(i-7)- 



In pQ < 1 limit [331, 



(43) 



(44) 



E'n-ip)^P 



1-27 



Therefore, the integral can be cast into the form 



Q27-lr2(3 26-127 



(pQ) 



27-1 



dh — integral 



^ +200 



dj 



PaPb ^ TT— loo 



2x(i) 



X(7) 



1 1 

53P 1 - 7, 77, 



r(7) 



2' 27 r(l-7) \p^q^ 



7-1 



(i7 



(1 1 1 

53P 1 - 7, o 



pIpIj\-,oo 27ri 2x(5) -x(7) 

r2 (I - 7) 26-127 ^4Q2y^-l ^p2g2 X7 



r(7) 



2' 2/ r (1-7) 



(45) 



(46) 



According to the residue theorem, the contour integral equals to the sum of all residues of poles enclosed by 
the contour. In the p^gi <C 1 limit, the first poles to the left and right of the vertical contour are the dominant 
contribution. The positions of these two poles are determined by the equation 2% (i) — x (7) = 0, and they are 



76 



(2)* 



0.18 and 1 — 79^'* = 0.82. For the first integral(see the left graph of Fig. (2]), we should close the contour to 



9 



the left to — oo and compute the residue of the pole at 7 = 0.18; while for the second integral(see the right graph of 
Fig. [5]), we should close the contour to the right to 00 and compute the residue at 7 = 0.82. Eventually, up to leading 
order precision, one reaches the result 

^Y,Q,q^,qi,Q - TT W + la + Qb ) riy^Q 

^^'^ iQ + qa+ qb) ■ 



1 ^^2^2x1-7^* 



"^PaQ-Pbq- \ 4 



93P ( 










X' 


bo'*) 





r (1 - 7^''*) (i + 7^'^*) 



r(i-7r) r(7r) r2(f-7r) 



^,2 



exp [2 (ckp — \)Y] cxp 



^e{p^,qa)e{ph,qb) ■ 



(7aC (3) CfYY 



(47) 



In this leading order calculation of 1 2 amplitude, we push the rapidity integration to its upper limit Y which leaves 
infinitesimal rapidity for the upper pomeron. As a result, the anomalous dimensions is determined by the dynamical 
pole introduced by the rapidity integration. This result agrees with eq.(44) in ref.j^. 



"2 (y, xoi , q) 



ttxqi exp [2 (ap — 1) y] 
8 (7aC (3) CpYf 



+00 



dv 
2^ 



Trexp [2 (ap - l)Y] [q'xf,^ 



-00 41n2-xH 



8(7aC (3)Cpy)' 



(48) 
(49) 



where Vi, is defined by 



and 



J-00 



d^X2 

XQIX12 



Jo -zqxoi e {xi2,q) e ixo2,q) , 



(50) 



(2)*' , 

7o = 2 + 



(2)* 



0.82. 



(51) 



In reaching the final result, we have assumed that has no pole before 7q . The proof of this assumption is 
explicitly provided in appendix ID] 



B. n 7^ part 

n ^ part of the amplitude is angular-dependent, and it vanishes after averaging over the angle of initial dipole 
orientation. This part of amplitude comes from the subdominant trajectories of the upper pomeron. As discussed in 
Appendix [B] we should restrict ourselves to the forward scattering case in which (5 = 0, since we are unable to prove 

that the whole discussion can be generalized to arbitrary momenta. According to the discussion above, one can easily 

(2) 

obtain the expression for Uy for nonzero value of n in forward scattering. 



''Y,Q=0,q- 



\n\>Q 



dh 



bn,u \P\ 



n/0 



pIpI 2x(i) -x(n,7)^ 



7-1 r3 



{h)T{l-2h) 
2'2)\ql) T{2h) r3(l-ft) 



1 PaPb 

47r2 



exp [2 (ap — l)Y] cxp 



eiPa^qa)eiPb,qb) 



28qCfC(3)V 



(7aC (3) CfYY 



(52) 
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FIG. 3: Illustration of the contours in complex 7 plane for the computation of ny when n 7^ 0. 



Changing / dh into J and using the identity gap (l ~ 51 5) — 93P (l ~ 5) which we proved in Appendix 



IfI along with the detailed derivation in Appendix [El one can reach 



''Y,Q=0,q- 



\n\>0 



n2 „2 



c?7 



(7+¥) 



x*l"l (p) + i?-"^" (p)] 53P (1 - K \, 

exp [2 (ap — 1) y] exp 

47r2 g2 



7-1 



^(ln''(';aPa)+ln^fePb)) 
28qCfC(3)V 



(7aC (3) Cpyr 



In pQ < 1 limit [331, 



iJ^^'' (p) = Q*'^-"/'q-+"/22-6- 



(53) 



(54) 



PQ 
4 



(-1) 



Q^^-^r(|-7 + f)r(|-7-§)2^-^^^ /QV 



(pQ)'"''. (55) 



r(i + 7 + f)r(i + 7-f) \QJ VP. 

where and parts give the angular dependence, and vanish after averaging. Thus, Eq''^ (p) + Eq "'"^ (p) — 

^1-27 Similarly, one has to examine the pole structure of the integrand of the 7 integral. Assuming 



that 53P (1 — h, i, i) does not contribute any new singularity in the domain ( — + 1 ) (see Appendix IF)) , noting 



that X ill IT') is analytic in the strip — < Rei' < + 1, and defining j^^'* as the solution to the equations 
2x (^) — X ("•, 7) =0 in the domain ( — 5 ) (see the left graph of Figl3]), one reaches 



I- 2 2^l-')'^^^'* 
oc e (p^, Qa) e (pft, gfc) 

Pal-PbQ- 



(2), 



exp [2 (ap — 1) Y] exp 



^(ln^(gaPa)+ln^(gi,Pl,)) 
28QCi?C(3)i' 



(7aC (3) CpF)^ 



(56) 



f 21 * 

by closing the contour to the left to —00 and collecting the residue at 7 = 7„ ' (other residues are suppressed by 
factors of p^q^). It seems that the dynamical pole introduced by the rapidity integrations always comes before other 
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FIG. 4: Dipole triplet density riy in an onium state with initial size p = pg — and three children dipoles having transverse 
separation p^, pj, and p^. The left graph shows the rapidity structure of ns, and the right graphs indicates the graphical 
representation of two successive triple pomeron vertices. 



/ 2 2 

(X -^^ — — still holds even for nonzero 



PaQ-PbQ- 



poles. Here we conjecture that, although we can not prove, ?^yQ=o'^'' 

momenta of Q by picking up poles at 7„ or 1 — 7„ depending on how the contour is closed(as shown in both 
graphs of Fig. [3]). 

(2) 

Comparing this result to the result in the case of n = 0, one can spot that ny with nonzero n is suppressed by 

factors of p^q^ in the p^q^ ^ 1 limit since 1 — 7!^''* ~ 1 + — e„, where e„ is a small positive number. Therefore, 
the sub-dominant trajectories can be neglected in the case of forward scattering computation. It is our conjecture 
that they can also be neglected in the non-forward scattering case. This result is new. 

To summarize, we have shown that the dipole pair density in momentum space scales as (p^f?^) " with respect 
to the onium size in p^q^ <C 1 limit. This indicates that the anomalous dimension of dipole pair density is equal 
to 0.18. Similar result can also be found in ref. ([Il])- Qualitatively, this result is easy to understand according 
to the discussion in Section |IT1 The anomalous dimension is determined by the upper pomeron which connects the 
triple pomeron vertex to the initial onium as seen in fig. [TJ The rapidity length of the upper pomeron is infinitesimal 
after rapidity integration. The anomalous dimension should approach zero in small rapidity limit or DGLAP limit 
{p^q^ <C 1 implies the coUinear ordering and the evolution in virtuality). In the case of onium-onium scattering, the 
double pomeron exchange amplitude between two onia with sizes p and p' then scales as 

1^(2) f / -xn (2 2\0-82 / 2 /9\0.82 f n CXp [(cKp - 1) F] 1 

t (7aC (3) CfF)^/^ J 



IV. THE 1 ^ 3 AMPLITUDE AND TRIPLE POMERON EXCHANGE. 



As shown in Fig. 21 the 1^3 amplitude is defined as the dipole triplet density in an onium state. It is easy to 
generalize the 1 =J> 2 amplitude to the 1 =J> 3 amplitude by adding one more triple-pomeron vertex. This amplitude 
corresponds to the graph in which one ancestor pomeron splits into two pomersons and one of the two descendent 
pomerons again splits into another two pomerons. As discussed in ref. [2^, the amplitude of this graph is different 
from the 1 3 amplitude introduced in ref. [1^ which involves a non-local 1 3 pomeron vertex(see also ref. (37|). 
As seen in ref. [28| . an initial pomeron can split into p {p > 1) pomerons simultaneously^ia a non-local 1 p pomeron 
vertex. The use of the non-local 1 p pomeron vertex {p > 1) is still under debate [29|. It seems more natural that 
the theory only requires one type of vertex (triple pomeron vertex) instead of infinite number of different vertices. 



Therefore, the 1 => 3 ampUtude reads, 

{PoPll PaoPai^ PboPbi^ PcoPci) 

dhdhadhbdhcdhd J d"^ p^d"^ p^d^ p^d^ p^d^ 



pIpIpI 

X 



I 



2m (J — (J {n, v) (J — uj (c) — ui {d) u> — ui {a) — (J (6) — a; (c) 

xE'^''^ip,SPu)E''-'^^ {PcoOPc^e) E^^^^ (PaoaPa.a) E^^^^" {p,,0P,J 

E^~^* {P28P,s) E"'-'"' {P2,P4,) e'"^'* {p,eP40) 



^J^2^^_P3^_P±T:^h,h*(„ „ \ T^hdJiA ( „ „ \-^h,,h,* 
\P23P34P42f 

IP56P57P67I 

Following the same procedure as in the last section, first of all, we define 

—h,l-ha,hb _ pl-h.h^.l-hc 



/ 



dy^d'p^dy ^,-,^ ^^^^^^^^ ^^^^^^^^ ^..rn.* (^^^^^^^ 



and 



-"■2,7,a,/3 — -"■2,7,a,/3 

\P56P57PeT\ 



I 



Using SL (2, C) transformation, one can easily reach, 

^?lA7f/''-'= = 9sp{l-h,h,,l-h.){p-r'Pes^''pJ,'^'') , 



where 



and 



i' = -h+1 - hd + he, 

f = -l + hd + hc + h, 
k' = -he + I - hd + h, 

i" = -hd + ha + h, 
j" = -ha + h + hd, 
k" = -hb + hd + ha- 



Moreover, in momentum space, we have 



/dhdhadhbdhcdhd [ 
J 2m uj — uj {n,v 



^) uj — uj {c) — uj {d) oj — oj {a) — ui (b) — lo (c) 

\pPaPbPc 



pIpIpI {2n^t 
X53P (1 -hd,l- ha, 1 - hb) gsp (1 - h, hd, 1 - he) 

X j Pad^ Ppd^ P-^d^ Ped^ Ps e^P {iQPs + kcPe + iQaPa + mPp) 
X {p-,l^' p-el^^' PS^^'') {T,l'^'p-e^'-~''p^^'-'') 

X {p-.r%T.^''' P-.I.-"''') (k^'p^^^"p:;^") 
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Hereafter in the computation of ny\ we will restrict ourselves to the limit of p^Q^ <^ 1 p^Qa ^ 1, ny, p^q^ ^ 1 and 
p^q^ ^ 1, where Q^, g^, q^ and q^ are of the same order. After changing of variables, 

w = \ {Pa+ Pp) + \ iPs + Pe): 

V ^ Pa- Pf3, 

W = Pa- p^+ Pf}- p^, 

X ^ Pe- Ps, 

V = Pe- P'y + Ps- Pi (70) 



the last three lines of Eq. ([69)1 can be written as 



J = —5^^HQ + qa + qb + qc) 



X j d?vd?w exp 
X /d^..^,exp 



// W — V 



.-1+/ f 



-i+fe" 



X [a.h.] 



x + y 



-l+i 



X [a./i. 



(71) 



where a.h. stands for the anti-holomorphic part of the square-bracket term, qab = qa — qb, qcQ ~ qc — Q and 
q'_ = I'^+'ii'-Vo-Q _ Following the same philosophy that wc have employed in the last section, we can expand 
exp \iq'_w — iqLy] into Taylor series (l + iq'_'w + • • • ) (l — iq'-y + ••■), and calculate the integrals order by order. 
One can easily see that all higher terms yield the same anomalous dimension as the first term does(See appendix lBjl . 
Thus, hereafter, we keep only the first term of the expansion 



Jo = ^S'-^\Q + qa + qb + qc) 

X (fvcfwexp i{qa-qb)T: 



, -i+i" / , \ -i+fe' 
-l+'i" I w — V \ / w + V 



X [a.h. 



X / (fx(fyexp i{qc-Q)7: 



-l+k' / , X 



X \a.h. 



(72) 



and then perform the integrals of coordinates 



TT 

Jo = -t^'-^HQ + Qa + qb + qc) 



qab 



2 \ +^ +'=-2 r (fc") r (j") r (^1 - J - fc j p (^" + f + k" - 1) 



—II —II 



r(fc" + j") 



r(i-r)r(i-r r 2-r-r-r 



2 \ '''+-''+'^'+* ^^^'^^^ r (A;') r (z') r ( 1 — / — 
qcQ 



r(z' + / + fc'-i) 



T{k' + 



r 1-? r 1-F r 2-?-J-I 



(73) 



In addition, at leading order in rapidity, one should pick the —^ 



iu—u{a)—u{b)—u{c) 



pole in / integral, and use saddle 



point approximation to evaluate / dhMK integrals which eventually fixes h. = h, ^ h. = I, = ^ ^ 0. 

Moreover, we assume n = = in the following discussion. (In Uy calculation, we have shown that higher n is 
suppressed. Here it is our assumption that higher n and would also be suppressed. Nevertheless, we have been 
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unable to find a proof for this point since the calculation becomes very lengthy.) Thus, 



n 



^(5^2' {Q + qa + qb + Qc) 



dhdhri 



3x1 



X (7) 2x 



1 1 



bo,u \p\ „o,i/ / \ (1 ^ ^ \ I 1 

4 V''"'r2(7,)r(i-7,)2i-4'^^ 



x47r 



x47r'^ 



ff2 
yafc 



7-7d 



1 e{Pa:qa)eiPh,qb)eiPc,qc) 



r2(i-7,)r(i + 7,) 

rM-7d + 1 + 7) r (I - 7d - 7) r (27, - i) 
r^Td + 5 - 7) r (i + 7d + 7) r (2 - 27,) 

exp [3 (ap — 1) y] exp 



8773 



Pa<labPbqabPcqcQ 



(74) 



Unfortunately, we are unable to to evaluate / dhd since is not necessarily a large parameter. This difficulty 

J '^"Q 

originates from the fact that all the relevant length scales have been integrated out. In this case we have to take into 
account all the singularities according to the residue theorem, or we can just choose the vertical contour from i — ioo 
to 5 + ioo and define the integral as a function of 7. Let us define 



dhd — integral 



did 1 



1 1 



.93P ( l"7d>2'2 '^^^ ( l"7,7d 



X (47r 



2^2 { qab_ 



r^(7.)rM-7. + ^ + 7)r(|-7.-7) 
r3(i-7,)r2 (7, + i-7)r(i + 7^ + 7) 



167r4^/ (7) 

IcQ 



(75) 
(76) 



Here one can not use the residue theorem since the ratio is not necessarily large or small. Nevertheless, it is 

straightforward to see that / (7) should be analytic in the strip domain < Re 7 < 1 when 7^ is integrated from 
i — ioo to i + ioo. Thus, 



3X(^)-X(7)^^^ (P)/(7)(^ 



7-1 



^ ^e{p^,qa)e{pi„qb)e{p^,qc) 



exp [3 {ap — 1) Y] exp 



^(ln^(gaPa)+ln^(gi,Pb)+'n^gcPj) 
28QCi?C(3)i' 



One can cast the final integral into 



(7aC(3)Cpr) 



9/2 



(77) 



,(3)P.p..P..P. ^ .(2).o + „ +„ /■d7 fii) 24" ^(i-^)^(^ + ^) |,iF°--r..i(' 4 ~ 

V,g„.,... -S (Q + ^« + '^^ + '^^)y^3x(i)-x(7)^ r(7)r(i-7) '^""^ ^"H^ 



7-1 



le(p„,ga)e(pfc,gb)e(pe,ge) 

X 



exp [3 (ap — 1 ) F] exp 



^(ln^(gaPa)+ln^(gi>Pt,) + lll^(gcPe)) 

28QCFC(3)r 



Pa1cQPb<lcQPcQcC 



(7aC (3) CpF) 



9/2 



(78) 



v Q '~ r { -—7)2 

In the end, following the same procedure as in the last section, we write Eq^ (p) o:^ p^~'^^ v'^['^^ ) — 



{pQ) 



27-1 



in pQ ^ 1 limit [33j, close the contour to the left for the first term, and close the contour to the right for the second. 
Therefore, according to the residue theorem, the contour integral equals to the sum of all residues enclosed by the 
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FIG. 5: Assuming 5y is infinitesimal, the anomalous dimensions of 1 => fc amplitude ny are illustrated in this figure. 



contour. The dominant contribution is 



■7o 



exp [3 (ap — 1) Y] cxp 



^{ln^{qaPj+ln^{qbPt)+ln^{q.p^)) 
28qCfC(3)>' 



Pa(lcQPb<lcQPc<lcQ [laC, (3) CpYf^ 

X(5^^^ (Q + + + qc) e (p^, 9a) e (p^, 9b) e (p^, q^) 



(79) 



where 7o^''* = 0.12 and 1 — 7o'^''* = 0.88 are the two solutions of the equation 3x {\) — X {l) — i'^i the domain [0, 1]. 

The anomalous dimension 7q = 0.12 found here is new. 

In the case of onium-onium scattering, the triple pomeron exchange amplitude between two onia with sizes p and 
p' then scales as 



(p,p',g,r)oc(9V)'' '' (9V^) 



0.88 ] 2 \{^P ~ 1) ^] 

(7aC (3) CpYf'^ 



a 



(80) 



GENERALIZATION TO 1 ^ it AMPLITUDE 



Based on the above computation, one can easily conjecture that the anomalous dimensions oil ^ k amplitudes 
should be governed by the equation fcx(^) — x{l) = 0, which comes from the dynamical pole introduced by the 
rapidity integration. Another way of stating this conjecture is that 1 ^ fc amplitudes tend to have constant energy 
dependence from to Y'(see Fig. [5]) as a result of kx (■^) = xil)- Namely, the total BFKL intercepts of multiple 
pomeron exchanges are independent of rapidity y. Supposing 7q'^''* is the solution to this equation between and i, 
then the 1 => /c amplitude should scale as 



(fc)(p,Pi---Pfc) 

I b-tr / \ KJ^ 



(p q 



-To 



exp [k (ap — 1) F] exp 



7T[ln''ipq)+ln^{p^qi)+-+ln''{p^qk)] 



(Piq) ■ ■ ■ (Pkq) (7aC (3) CpF)^''/^ 

X(5^^' {q + qi-\ hgfc)e(pi,9i) • ■■e{pk,qk) ■ 



(81) 
(82) 



For reader's convenience, we list the first five 7q in the following: 7q 



1 (2)* 
2' '0 



0.18, if* 



0.12, 7^,'*^* = 0.090, 

7g ' = 0.072. This result fits the naive expectation that the anomalous dimension should approach zero in small 

rapidity limit or in DGLAP limit. Furthermore, when k becomes large enough, we obtain 7o''''* ~ 0. This indicates 
that the anomalous dimensions of large order pomeron loops are dominated by the DGLAP evolution. We should 
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remind the reader that this conjecture can only be used in the region 

?i 9i ■ ■ ■ I 9fe are all of the same order, and so are all the dipole sizes p^p^ - ■ ■ ^ pj.. In the case of onium-onium scattering, 
we obtain that the corresponding k-pomeron exchange amplitude should scale as 

p(fc) / I v\ / 2 2\l-7;'>* / 2 /2\l-7j,'** / 2 exp [(ap - 1) r] 1 

^^()(p,P,,,r)oc(,p) I" (7«c(3)c.yf/^l ' ^''^ 

where we restrict ourselves to the p^q^ ^ 1 and p'^g^ ^ 1 limit. As a result, different order of pomcron-loop 
amplitudes have different dipole size dependence, namely, they belong to different universality classes. Thus, it seems 
that re-summation of leading order amplitudes is insufhcient and impossible to reach unitarity. The argument is 
straightforward: supposing that unitarity is achieved by resummation of all loop amplitudes by means of delicate 
balances for fixed p, it seems impossible to obtain unitarity again when p is changed to another value since the 
delicate balance is then broken. 

In summary, according to the conjectured formula, we have found infinite number of new anomalous dimensions 
between ^ and 0. Comparing the multiple porneron exchanges to single pomeron exchange, it seems that the evolu- 
tion gets pushed towards DGLAP evolution 34| , and the anomalous dimension discretely approaches the anomalous 
dimension of DGLAP evolution. Because of the lack of the rapidity space for upper and lower pomerons, the evo- 
lutions in pomeron loops manage to balance themselves in between the BFKL evolution and the DGLAP evolution. 
This explains why the new anomalous dimensions are distributed in between i and 0. Furthermore, it is intuitive to 
notice that the saturation anomalous dimension found in ref. [ssj ] is 7^ = 0.37 as a solution to Kovchegov equation 
in geometric scaling region. Kovchegov equation essentially resums multiple pomeron exchanges (fan diagrams) and 
yields an anomalous dimension between i and as a result of re-summation of multiple pomeron exchanges. The 

explicit connection between the 7j, and 7q , however, is still unknown and remains as an open question. 



VI. COMPARISON WITH REGGEON FIELD THEORY CALCULUS 

Reggeon field theory (RFT) calculus(for a review, see ref. [1^, EO, lU, ) , similar as Feynman rules, provides definite 
intercepts {a[t) — 1 = a (0) -I- a't) and propagators for reggeons (including pomerons) in QCD phenomenology. In 
RFT, the pomeron anomalous dimension 7 = i is a universal (conserved) quantity and it corresponds to the pomeron 
intercept ap ~ 1 — Is£Z£iin2_ leading order of 1 ^ 2 amplitude, RFT has a genuine triple pomeron vertex which 
connects 3 pomerons with anomalous dimension 7 = ^- This approach certainly is justified in some situations such 
as large diffractive mass scattering, where the large diffractive mass is large enough to fix the anomalous dimension 
of upper pomeron at \ by saddle point approximation(e.g., see ref. Q). On the other hand, as far as the unitarity 
problem is concerned, one has to integrate the intermediate rapidity to the upper limit Y , in which case the anomalous 
dimension of the upper pomeron in the vertex is no longer fixed at i. 

Following this philosophy and the essence of the RFT, in order to compare with the calculation we have finished 

f2l 

above, let us re-consider ny in forward scattering (Q ~ 0,qa — — <z& — q) by assuming that y — y is always large 
enough to justify the saddle point approximation (In fact, this assumption breaks down when one integrates rapidity 
to the upper limit Y .) First of all, let us begin with the expression with the rapidity integration. 



,(2)P:Pa 



pIpI 







bn. 



P\ Tpn,ii 



dy / dh exp [u (7) {Y - y)] 



7-1 



PaPb 



exp [2 (ap — 1) y] exp 



■47r2 



e{Pa,(l)e{Pb,-q) ■ 



r3 (7)r(i-27) 
r(27)r3(i-7) 

2SaCFC{3)y 



{7aC (3) CpyY 



Dropping all the higher n, and changing / dh into 



one can simplify riy^ and get. 



n. 



d7 fq^p^y-' . , ..^ 4r(7) 



exp [2 (ap — I) y] exp 



47r2 q'^p^p^ 



e{Pa,q)e{pb,-q) ■ 



1-7- - 
^'2' 2 



^(ln"(<?pj+ln"(gp,)) 
28qCfC(3)2/ 



(84) 



(7aC(3)Cpy)^ 



(85) 
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If one uses saddle point approximation to evaluate the J ^ integral before the rapidity integral, one reaches, 



"'y,Y,Q,g-,0 



Stt^ qPaQPb \2 2 2J 
exp [2 {ap -l)y+ (ap - 1) (F - y)] exp 



7i-(ln''(gp„)+ln''(gpt,)) 7r(ln''(i;p)) 

28aCFC(3)y 28aCFCi3)(Y-y) 



(7«C (3) Cpy)' [7aC (3) Cp (F - y)f' 



(86) 



Indeed, we now obtain the genuine triple pomeron vertex which connects 3 pomerons with anomalous dimension 
7 = i. Finishing the rapidity integration yields. 



e\ 93P (Ll^l ) e(Pa:'?)e(Pb,-9) 

QPaQPb V2 2 2.' 



-exp [2 (ap — 1) Y] exp 



7r(ln"(gpJ+ln"(gpt,)) 



28qCfC(3)>' 



(87) 



(88) 



Vl4C(3)ln2 (7aC (3) Cpy)' 

This result agrees with eq.(56) in ref.Q. Integrating the triple pomeron vertex found in eq.(56) in ref.[5'], one gets 



yn2 {Y,y,XQi,q) 



Y 



dy- 



yCfVoXoi exp [2 (ap -l)y + (ap - 1) (y - y)] 



8<zV7«C(3)Cp(y-y) (7aC (3) Cpyf 

nVo exp [2 (ap - 1) Y] xqi 



16 V14C (3) In 2 (7aC (3) Cpy)"* 9 



(89) 
(90) 



It clearly differs from Ea. (|T7|) which is proportional to 

J2) 



-53P (0-82,^,^) as a result of the anomalous 
dimension being 1 — 7q^^'. The origin of this difference comes from the rapidity integral which dynamically changes 



PaQ-Pbl- 



the anomalous dimension from h to 



2 To • 

The zero transverse dimension toy models [3 |4^, \^ [4^, |4^, which catch the essence of reggeon field theory 
calculus, also have universal rules for the pomeron intercepts and propagators, as well as the triple pomeron (or 
reggeon) vertices. The toy model does not contain the anomalous dimension or transverse dimensions. In some sense, 
it over-simplifies the problem and fails to catch new features of the QCD pomeron that we found above. 

On the other hand, the QCD dipole model, as one can easily spot from the calculation in Sections |TT1 IIIII and 
rvl contains the distinct feature of the non-existence of universal intercepts and anomalous dimensions for QCD 
pomerons in the unitarity calculation. In leading order calculation of 1 fc amplitude, we always push the rapidity 
integral to the upper limit Y which leaves infinitesimal rapidity for the upper pomeron (see Fig. [S]). As a result, the 
anomalous dimensions of this pomeron are then fixed by the dynamical pole introduced by the rapidity integration. 
The anomalous dimensions are no longer universal and they vary from graph to graph according to the detail structure 
of the graph and dynamics. Nevertheless, the other side of the coin is that we now have the constant energy dependence 
(coefRcients of the rapidity) from to F while this certainly is not true in RFT. 



VII. CONCLUSION 



We explicitly calculate the anomalous dimensions of the single, double pomeron and triple pomeron exchange 
amplitudes in QCD dipole model in the leading logarithmic approximation. Other than the usual pomeron anomalous 

dimension ^, we find 7q =0.18 and 7q =0.12 for double and triple pomeron exchange amplitude, respectively. 
Based on the calculation, the general formula 7o'°''* = x^^ {^X (5)) which governs anomalous dimensions of 1 /c 
amplitude is conjectured in the region In^ (qp) ^ i4aCFC(3) y a^d p^q^ <C 1, where p stands for all the dipole sizes 
and q represents all the momenta scales. 

(2) 

Furthermore, the calculation of forward scattering ny in the ri 7^ case shows that contributions of sub-dominant 
trajectories of the upper pomeron are suppressed by powers of p^g^ in p^g^ <C 1 hmit. It is our conjecture that sub- 
dominant trajectories can be neglect in the non-forward scattering case. We utilize this conjecture in the calculation 

(3) 

of ny and get rid of the contribution from nonzero value of Ud- 
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In addition, different pomeron loop amplitudes (leading order) belong to different universality class as a result of 
different anomalous dimensions. It seems that re-summation of these amplitudes is insufficient and impossible to 
reach unitarity. One may have to take higher order contributions into account. 

Last but not least, in comparison with the reggeon field theory, one finds that there are two differences between 
this computation and the reggeon field theory although the BFKL characteristic function x il) is universal. The first 
one is that the anomalous dimensions are no longer a constant and they vary according to their positions in the graph 
while 1 = \ everywhere in RFT in order to have a fixed pomeron intercept (for small and fixed t value); the second 
difference is that the QCD dipole model tends to have a constant energy dependence (coefficients of the rapidity) from 
to y while this certainly is not true in RFT. Namely, the total BFKL intercepts of multiple pomeron exchanges are 
independent of rapidity y in QCD dipole model. The constancy of the energy dependence is equivalent to the general 
formula 7q'^''* — (^X (5)) can easily explain its physical meaning. 
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APPENDIX A: THE EVALUATION OF INTEGRALS. 

In this and the next appendices, we evaluate the following integral, 

where a.h. stands for the antiholomorphic part of the square-bracket term. First of all, let us perform the Taylor 
expansion of the exp (— iQ-j) , and keep only the first term 1. We will discuss the results for higher terms in the next 
appendix. Thus, the above integral becomes. 



1. The d^w integral 

Furthermore, one can decouple Ea. (|A2[) into two 2-dim integrals in which the d^w integral can be reduced to the 
following integral: 

/ = / ^.-1+^-1+^ (1 - (1 - zf-' , (A3) 

where i — i and j — j are integers. The solution to this integral can be found in Dotsenko and Fateev's paperjs^l in 
statistical physics. In the following, we carry out the detailed calculation (5ll| in complex plane. 

The first step is to change z = x + iy,'z = x — iy, and then perform a Wick rotation j/ — s- i (1 — 2ie) y, where the e, 
which is an infinitesimal positive number, makes sure that the singularities are not touched. We obtain, 

/+00 r+oc _ _ 

dx / dy{x — y + 2ieyy {x + y — 2ieyy^ {1 — x + y — 2ieyy {1 — x — y + 2ieyy . (A4) 
-00 J —OQ 

Next, one can change the variables into AT^ — x + y and AT- ~ x — y and cast the integral into the form. 



-l+i 



-1 + fc 



\a.h. 



(Al) 



-l+i 



-1+fc 



\a.h. 



(A2) 
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Jf+ < 

1 

^AAA^ Ay\_AAAy 



0<X+<1 



1 

AAAAAy 



X+ > 1 



^AAAJ• ^vww 



1 



+0O X- 



FIG. 6: Integral contour of X- and branch cuts in three different domains of X+, where the straight hne stands for the contour 
and the photon hne represents the branch cuts. 

0<X+<1 ^ 
C 



+00 x_ 



FIG. 7: The contour C. 



/ = / dX+ dX^ [X_ + ^e {X+ - X^)]'-^ [X+ - le {X+ - 

2 J— oo J —oo 

x[l-X^ -ie {X+ - X^)Y^^ [1-X++ ie {X+ - X^)Y^^ . (A5) 

The integral over X+ can be decomposed into 3 pieces with respect to three integration domains (see Fig. [5]). Only 
in the second domain does the integral give non-trivial contribution, since the contour of X^ can not be deformed to 
a single point only in that case. Therefore, we can reach the factorized integrals, 

/ = /' dX+xl-' (1 - X+y-' I ^X!_-' (1 - X^y-' , (A6) 

where contour C(see Fig. [7]) is a contour which encloses the branch cut [1, oo]. It starts from ie + oo, goes to ie + 1, 
then crosses the real axis to — ie + 1, in the end, goes to — ie + oo and forms a close contour. It is straightforward to 
compute this contour integral and obtain the final result, 

/ = [\x+xl-' {1 - x+y-' 

Jo 

X [exp (-z^ (j - 1)) - exp {itt (j - 1))] / \-X!_-' (X_ - if^^ , (A7) 

Ji 2i 

^ r(»)r(7) Tji-t-j) ^ T{^)^{J) r(i-z-J) 

r(z + J) r(i-i)r(i-j) r(* + j) r(i-z)r(i-7) ^ ' 



20 



Therefore, the cfw integral yields, 



/ 



w — V \ I W + V 



la.h.] 



= , /(fc)r(j) r(i-fc-j) 



T{k + j) r(i-fc)r(i-i) 



2. The d^v integral 



The (fv integral can be cast into 



/ 



oo p2Tr 



/ dr dipr^^^ ^ exp [i (77 — ry) ly? + ig_r costp] 
Jo Jo 



ri+r}-2 



(A9) 
(AID) 
(All) 



where rj = i + j + k and 77 — 77 = — (n + ria + ) S ^- In reaching the above result, we have set the orientation of the 
'q- parallel to the real axis^(^g_ = 0^, and used the following two formulae: 



Jm (z) 



1 



27ri™ Jo 



cxp {iz cos ip + imif) . 



/ dxxi^T (ax) - — V2 ^ 2 ^ 2; 

Jo " ^ ^ ~ a>^+' r (i + f - f ) 



(A12) 
(A13) 



where m should be integers. 

Combining the two results, one gets, 



jh,hM^^^2--in+n.+n,) / J_y+^''+^''-^ T {h + K -_hb)r {h + h - F (l - 2h)_ T jh + K +_hb -I) 



'Q,q- 



r (2/1) r{i-h-ha + hb)r{i-h-hb + ha)r{2-h-ha-hb) 

(A14) 



APPENDIX B: OTHER TERMS IN THE TAYLOR EXPANSION 

In this part, we discuss cases of the higher terms of the Taylor series of exp (iQ^). 



1. The second term in the expansion 

The relevant integral of the second term in the expansion reads. 



^Q,q- d vd wexp 



{iq- ■ v) 



w — V \ I w + V 



la.h.] i — • 
^ ^ 2 



w — v\ / W + V 



I' 



Q 



d^vexp {iq- ■ v) ( i— ■ v] '^Tf "^Att 



r(fc)r(j) r(i-fc-j) 



(Bl) 
(B2) 



27- " '■■r(i + A: + j)r(i-fc)r(i-J)" 

where 7^'^^<» stands for the relevant integral of the I + 1 term in the Taylor expansion. Assume ■ v = 
i%vcos{(j) + ip), where (j) is the orientation of Q in 2-dim plane, (j) should also be considered as the angle between 
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~q- and Q. It is straightforward to compute the v integral and obtain 

Q 



Q'1- 2g_ 



r(i + fc + i)r(i-^)r(i-J) 



r(H-i + fc-i) r(i + i + fc-i) 

exp (-19) = — = — =^ - exp [i-^^ 



r(i 



k) 



r(i 



k) 



(B3) 



Since k — j = 2 {ha — h^) and eventually ha = h}, = \, jh,ha,ht„i y^jjjgjjgg jj^ ^j^g g^^j Moreover, it is easy to prove 
that integrals of all odd power of are proportional to k — j and they vanish as well. 



2. The third term in the expansion 

The relevant integral of the third term in the expansion reads, 



I 



d'^vd'^w exp (zg_ • v) 
d^v exp (iq 



w — V 



-i+j / , \ -i+fe' 



[a.h.]^\^^■w 



(B4) 



x47r- 



2 V 2 

r(fc)r(j) r(i-fc-7) 



{k-jf + {k + j) 



T{2 + k + j)T{l-k)r{l~j 
Dropping the {k — j)^ term and finishing the J dPv integral, one geti 

IKHM = ,^2..-(„+„.+n.)QV4\W^--^ T{k)T{j) T{l-k-D 



(B5) 



'Q,q- 



r(r?) 

r(l-r?) 



T{2 + k+j)T{l-k)T{l-3) 
exp (2i^) r (77 + 1) exp (-2i^) r (?7 - 1) 



(k + j) 



T{2-rj) 



r(-^) 



(B6) 



Thus, it is straightforward to examine the singularity structure of ny^Q g'^'^'' {it is proportional to j^'^'^'^'"^ which 
comes from the third term in the expansion) before the final 7 integral in the complex 7 plane when n = 0. 



(2)P,P„,P6 



OC 



2m TT* 2x(i) -X(7) 

dl \P\ rrO,u 



E^^^{p)g,P [1-^,1^ 



2^z2x(i)-x(7) "'"'^^ V "2' 27 7+, 
In addition, one can also show that the fifth term in the expansion 



Q,q- 

r(7 + i) 
r(-7) 



cos 20 



r(7 + 2) 
r(i-7) 



(B7) 



OC 



\P\ 



I 



2m 7r4 2x(i)-x(7) 
dl \p\ 



2m 2x(i)-x(7) 



4'^(p)53p(l-7,^,^ 



1 



(7+ I) (7+i) 



(B8) 



where fi [7, </>] are terms which are singular only at integer values of 7. Therefore, we can see the pattern and conclude 
that 7 = 0.18 and 0.82 are the only two singularities in the (0, 1) domain when n = 0. In the case when n 7^ 0, since 
things are much more complicated, we have been unable to reach similar conclusion. Therefore, wc have to restrict 
ourselves to the forward scattering case in which Q = 0. Nevertheless, contributions from sub-dominant trajectories 
would not play any role when one averages over the orientation of the initial onium, since those contributions are all 
angular dependent and vanish after averaging over angles. 
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APPENDIX C: SADDLE POINT EVALUATION OF THE POMERON 



In this appendix, we use saddle point approximation to evaluate the pomeron trajectory when Y is large. The 
relevant integral is 



27r2 



+ 00 

J J (fRex_j>{iqa-R) 



1+211/ „ 



\R-^\ \R+^\q 



exp [uj (a) Y] , 



-\-oo 



exp (zga • R) 2/9„ /" di^a 



d Ri o 11 01 / ^2;.^ exp 



2 In 2asY - 7^ (3) lyla^Y + 2iva In ■ 



2Pa 



(CI) 
(C2) 

, (C3) 



where as = ^"'^"^^ ■ In reaching the above result, we have used the idea of saddle point approximation and neglected 
all sub-dominant trajectories. In the region where the diffusion approximation In^ (qp) il£C^£C(3ly jg valid, this 
integral yields 



1 exp [{ap ~ 1) Y] exp 

^7" i7aC{3)CFYf' 



Trln (p„i?a) 
28qCfC(3)>' 



(C4) 



where we have set In 



2Pa 



Id Pa Mr PaV ~ (PaQa) since J (PR is dominated in the region where i? ~ l/^ai a-nd we 
have defined e (p„ = 27 / t^'^exp i ^^p^L^^p^, , ap - 1 = sin^^. 



APPENDIX D: EXPLICIT CALCULATION OF K 

According to Eg. ([50)1 and Mellin transform, one can easily obtain. 



(fx2 



1 



^ {I'^^oiy ^ ^7 = / -^0 ( ^'?a;oi ) e (2:12, g) e (xo2, g) ■ 



(Dl) 



Evaluating the above integral in g^Xgi <C 1 limit, we should close the contour to the right half plane. The leading 
order term tells us the position of the first pole of Vj in complex 7 plane. 
First of all, it is straight forward to discover that 



1 exp \iq ■ X (i — u)] 
e{x,q) = - du^^ ^^TTTT^^O 



Jo [u{l-u)\ 



1/2 



qx\J u{\ — u) 



where we have used the formulae 



(D2) 



du 



and 



(x^ + a^) 



Changing variable u — sin^ |, and combining 



J, (qx) x''+^ ^ a'^-^qt'K,.^ {qa) 
rdx — 



2\/^+i 



2'^r (m + 1) 



2sin7r/i ^ ^ 
If, (x) = exp Jf^ exp (-«^) 



(D3) 

(D4) 

(D5) 
(D6) 
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along with the integration identity 



r/2 



[z sin t] cos (x cos t) dt = i^Ji^ [v+\ Ji^ [U-] ■, 



where y± 



gives 



e{x,q) 



hm — 

Ai^o 2 smTT/i 



J_.(^.)J_. (e*)- J^(^>) J^(^^*) 



where g = ^e"^ with being the angle between ~q and ~x . Therefore, 

/ ^ Jo 1 ^^a^oi ) e (a;i2 , g) e (2:02 , g) 
7 X01X12 \2 J 

= '^Jo (^Y^) j dxi2dx2QdbbJo (bxoi) Jq (&X02) Jo (^2^12) 

X [^0 (£•12) Jo (£'^2) + Jo (£'12) ^0 (£'12)] [^0 (£'20) Jo (£'20) + Jo (£'20) Yo (£'20)] > 



where the identities!: 



and 



d'^X2 = J {X21,X20) dx2idx2 



- I dbbJa {bxoi) Jq {bxo2) Jo {bxi2) = 



4a;2iX2o 



{X21 + X20) - xIq - {X21 - X20) 



1 



(X21+X20) - xIq - {X21 - X20y 



(D7) 

(D8) 
(D9) 



(DIO) 
(Dll) 

(D12) 



have been used. It is easy to estimate the above integral and obtain f -^-^^Jq (hqxoi) e (xi2, q) e (xo2, q) In -jt-, 
which indicates that the first pole of occurs at 7 = 1. 



APPENDIX E: DETAILED CALCULATION FOR HIGHER POMERON TRAJECTORIES. 



Here we consider the integrand of J d'y in the case of general value of ; 

2 iv"^ + n^/A) hh h 
Integrand = 2tt bn^Id *" 



4(i^2_^n4/4) 71324*"^ T{\n\/2-iv + l/2)T{\n\/2 + iv) 
\n\ /2-ivT {\n\ /2 + ii' + 1/2) V {\n\ /2 ~ iv) 

,,,.,-.,Y4^'^3(M^(l-2E) 



gV T{2h)T^{l-h) 



(El) 



(E2) 



When n > 0, the integrand yields (-i)"'"' Svr^ (^^^ r(i+M'!,^) ' ^"^^^ < 0, it gives ^'"ISTr^ (^^^ Jlz^^iL. 



Generally, it can be cast into zl"l87r'^ r(i+"i"i'' ^ ) ' 



^(1+^-7)- 



APPENDIX F: TRIPLE POMERON COEFFICIENT qsp 

In this appendix, following Korchemskv [3]| . we discuss the triple pomeron coefRcient g^p. 

72 „ ^2, 



53P ih,ha,hb) 



d^P2d P^d^PA h -h ( P2A 



P23P34P42I 
P34 



P2P4 



(1-P3)(1-P4); V(l-P3)(l-P4) 



ha / — \ ha 

P2A 
P2P4 



P34 



(Fl) 
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Changing the variables: p' = p and p' = p, 



P34 



P34 



/lb 



(F2) 
(F3) 



which means triple pomeron coefRcient g^p is an even function of n, Ua and Ub. Following the methods used in 
ref. [3l[ (here our h corresponds the h-^ in ref. [3l[), we find that 



(F4) 



where Jc [h, i, i) and Jc [h, i, i) can be written in terms of Meijer' G functions G|| and hypergeometric functions 
4F3 as follows 



1, i + /l, 1, I - ft.' 

1111 

2' 2' 2' 2 , 



= F (1 - h) Gli 1 



1111 

2 ' 2 ' 2 ' 2 

i-ft,0,ft-i 



(F5) 



J2 h 



1 1 

2' 2 



J3 h, 



1 1 



n^T^ (1 - ft) 



2'2; F(f-ft)F(|-ft) 



4F3 



i,i,l-ft,l-ft 
1, 1 - ft, I - ft 



1 ; 



(F6) 



J3 (ft, fta, ftb) = J2 (ft, ftb, fta) 



(F7) 



Ji f ft, — , — 
V 2' 2 



7r2 



F(l-ft)F(ft) 



1, I -ft, 1,1 

1 1 1 _ 1 _ 
2'2'2 2 ' 2 . 



7r2 



=^ ^=rG^^ 1 



F (1 - ft) F (ft) 



1 1 17^+ 1 17^' 
2 ' 2 ' 2 ^ 2 ' 2 

0,ft-i,0,0 , 



(F8) 



42 
44 



1, § - ft, 1, i + ft 



1 

2 ^' 2 

1 1 1 1 

2 ' 2 ' 2 ' 2 



^24 
lj-44 



1 1 1 1 

2' 2' 2' 2 

0,ft- i,0,i -ft 



(F9) 



J3 (ft, fta, ftb) = J2 (ft, ftb, fta) 



(FIO) 



We find different but equivalent final expressions of Jc (ft, fta, ftfc) and Jc (ft, fta, ftfc) as compared to those in ref.[3l|. ^ 
Numerically, one finds (5, ^, ^) = 7766.68 which agrees with the previous result in ref . [3l|, [l^l ■ Thus, we outline 
our evaluation of the g^p (ft, ^, ^) in the following. For instance, from ref. [sH] . one can get 



F (1 - ft) F (ft) Jo 



2^1 



dxx ^ {1 — x)" 



I I 

2' 2 
1 



X I 2F1 ( Y 



Using the identitvfs^ 



(Fll) 



2i^l 



1 1 

2' 2 

1 



X 2F1 



I 1 
2' 2 

1 



{l-x) 



1/2 



3^2 



111 

2' 2' 2 

1,1 



4(1 -.t) 



(F12) 



^ Due to some numerical subtleties, B.X. did not realize the equivalence between the results obtained above and the final expressions of 
Jc {h, ha, hi,) and Jc (h, ha, /»(,) in ref. [3l[. B.X. would like to thank Dr. L. Motyka for communications on this issue. 
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and the Mellin-Barnes representation of hypergeometric function 

1 /•+*°° ds r3(i + s)r(-^ 



3-^2 



111 

2' 2' 2 

1,1 



along with identities 



Jo 

r(2s) 



2771 r2 (1 + s) 



r(/x + i.) ' 

22S-1 



-r(s)r s + 



one obtains, 



^- 1 1 
■^i|'^'2'2 



7r2- 



r (1 - /i) r (h) 

'■+'°° ds r2(i + s)r(s + i-p)r(s + i-p)r(-s)r(7^-i-s) 



27ri 



According to the definition of the Meijer's G function 



r2 (1 + s) 



f^mn 

pg 



tti, . . . , ttp 
bi,...,bq 



n r(i-&i-s) n r(a, + s) 



j=n+l 



(F13) 

(F14) 
(F15) 



(F16) 



(F17) 



one can easily reach the expression wc found above. It is straightforward to put g^p (/i, ^, ^) into Mathematica and 
plot it for fixed integer values of n as function of 7. The plots for the first few integer values of n show that gsp (h, ^, |) 

has no pole in the domain ^— -ly-, ^ + 1^. 
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